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The recently discovered r-mod e instability ( Andersson 



1998a 



has significant implications on the rotational evolution of 
a newly-born neutron star. The r-modes are unstable due 
t o the Chandrasekhar - Friedman-Schutz (CFS) m echanism 
( |Chandrasekhar 1970^ |Friedman fc Schutz 1978| ). Two in- 
dependent computations by Andersson, Kokkotas & Schutz 
(1998) and Lindblom, Owen & Morsink (1998) find that the 
r-mode instability is responsible for slowing down a rapidly 
rotating, newly-born neutron star to rotation rates compa- 
rable to that of the initial period of the Crab pulsar (~19 
ms) . This is achieved by the emission of current-quadrupole 
gravitational waves, which reduce the angular momentum of 
the star. The instability is active for as long as its growth- 
time is shorter than the damping-time due to the viscosity 
of neutron-star matter. The r-mode instability explains why 
only slowly-rotating pulsars are associated with supernova 
remnants. The r-mode instability does not allow millisecond 
pulsars to be f ormed after an accretion -induced collapse of 
a white dwarf (Andersson et al. 1998a). It seems that mil- 
lisecond pulsars can only be formed by the accretion-induced 
spin-up of old, cold, neutron stars. Additionally, this insta- 
bility should be active in accreting neutron stars and so can 
limit their rotation provided that the stars are hotter than 
about 2 x 10 5 K and also in the rapidly spinning neutron 



ABSTRACT 

We present a mathematically rigorous proof that the r-mode spectrum of relativistic 
stars to the rotational lowest order has a continuous part. A rigorous definition of this 
spectrum is given in terms of the spectrum of a continuous linear operator. This study 
verifies earlier results by Kojima (1998) about the nature of the r-mode spectrum. 
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as a sum of spheroidal and toroidal components (or polar 
and axial components, in a different terminology). In a non- 
rotating star, the usual /, p and g modes of oscillation are 
purely spheroidal, characterized by the indices (l,m) of the 
spherical harmonic function Y; m . In a rotating star, modes 
that reduce to purely spheroidal modes in the non-rotating 
limit, also acquire toroidal components. Conversely, r-modes 
in a non- rotating star are purely toroidal modes with van- 
ishing frequencyrj 

In a rotating star, the displacement vector acquires 
spheroidal components and the frequency in the rotating 
frame, to first order in the rotational frequency Q of the 
star, becomes 



Friedman & Morsink 199S ) in rotating neutron stars, 



2mQ. 
1(1 + 1) 



(1) 



for a given (I, m) mode. An inertial observer, measures a 
frequency of 



mQ, 



(2) 



stars i n low mass X-ray binaries (LMXB)( Andersson et al 



1998b). Finally, while the initially rapidly rotating star spins 



down, an energy equivalent to roughly 1% of a solar mass is 
radiated in gravitational waves, making the proce ss an in- 
teresting source of detectable gravitational waves (Owen et 
al. 1998j ). 

Oscillations of stars are commonly described by the La- 



As mentioned earlier, when the star is set in slow rota- 
tion the axial (toroidal) modes are no longer degenerate, but 
instead the new family of r-modes emerges, which are hori- 
zontal displacements on equipotential surfaces. In this case 
the axial (toroidal) perturbation with spherical harmonic 
indices (£, m) induce polar (spheroidal) perturbations with 
harmonic indices (t ± l,m) and vice versa. 

The picture described above is purely Newtonian, the 
above calculation was performed using Newtonian slowly ro- 
tating stellar models and the power radiated in gravitational 
waves was estimated using the quadrapole formula. The two 



grangian displacement vector £, which describes the dis- 
placement ot a given fluid element due to the oscillation. 
Since £ is a vector on the (6, <j>) 2-sphere, it can be written 



In the rela tivistic case the picture is similar ( Thornc & Campo 



lattaro 1967P nevertheless, in this case there exist an additional 



family of quasiiiormal modes, the ones called "svacetim e or w- 
modes" ( Chandrasekhar & Ferrari 1991b; Kokkotas 1994). 
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approximations (Newtonian theory and slow rotation) that 
were used give only a qualitative picture while the quanti- 
tative results would change if at least general relativity is 
used. The assumption of slow rotation is a robust approxi- 
mation because the expansion parameter e = Q.yJ R 3 JM is 
usually very small and the fastest spinning known pulsar has 
e~0.3. 

The perturbation equations for slowly rotating rel ativis- 



tic sta rs were derived by Kojima (1992), see also (Chan- 
drasektar & Ferrari 1991a). Kojima (1993) calculated also 



the effect of slow rotation on /-modes. Andersson (1998) 
found the r-mode instability using the same set of equations, 
although his calculations overestimate the growth rate of the 
instability. 

An important difference between Newtonian and gen- 
eral relativistic calculations is the dragging of the inertial 
frames, which might produce significant changes in the fre- 
quency spectra. This is exactly the case that we a re studying 
in this article. To be more specific, Kojima ( 199S ) suggested 
that if one calculates the r-mode frequencies using general 
relativity to lowest order in Q the spectrum becomes contin- 
uous, in contrast to the calculations from Newtonian theory 
where the spectrum is discrete and the frequencies are given 
by the formula (Q) . In this article we prove in a mathemati- 
cally rigorous way that Kojima 's suggestion for the existence 
of a continuous part within the spectrum is true. 

Continuous spectra have been found in many c ases in 
the pas t in the study of differentially rotating fluids (ffohut? 
& Verdaguer 1983[ [Verdaguer 1983;; iBalbinski 1984al |Bal- 



binski ; 9841 ) . The continuous spectrum in these cases were 



again seen for the r-modes together with a wealth of inter- 
esting features such as: the passage of low-order r-modes 
from the discrete spectrum into the continuous one as the 
differential rotation increases; and the presence of low order 
discrete p-modes in the middle of the continuous spectrum i n 
the more rapidly rotating disks (Schutz & Verdaguer 1983). 



The stars under consideration here have no differential rota- 
tion and the existence of a continuous part of the spectrum 
is attributed to the dragging of the inertial frames due to 
general relativity. 



2 PERTURBATION EQUATIONS 

Since our calculatio n s will be based on the equations of Ko- 
jima (Kojima 1998; Kojima 1997) which are presented in 
detail there, here we are going only briefly to describe the 
perturbation equations. 

We assume that the star is uniformly rotating with an- 
gular velocity il ~ O(e) where e, as stated earlier, is small 
compared to unity. The metric is given by: 



ds 2 — — e v dt 2 + e x dr 2 + r 2 (d6 2 + sin 2 

n 2-2 

— 2wr sm 



(3) 



where u) ~ O(e) describes the dragging of the inertial frame. 
If we include the effects of rotation only to order e the config- 
ura tion is still s pherical, because the deformation is of order 
e 2 (Hartle 1967). The star then is described by the standard 



Tolm an-Oppcnheimer-V olkov (TOV) equations (cf Chapter 
23.5 (Misner et al. 1973)) plus an equation for u> 



where we have defined 
cc = f2 — LO 

a prime denotes derivative with respect to r, and 

In the vacuum outside the star w can be written 
o 2J 

W = 52 — . 



(5) 



(6) 



(7) 



where J is the angular momentum of the star. The func- 
tion w, both inside and outside the star is a function of r 
only and continuity of zu at the boundary (surface of the 
star, r = R) requires that zu' R = 6JR~ 4 . Additionally, zu is 
monotonically increasing function of r limited to 



Wo < ZU < il, 



(8) 



where zuq is the value at the center. 

For the study of the perturbations of slowly rotating rel- 
ativistic stars we should expand all perturbation functions 
in spherical harmonics and additionally assume a harmonic 
dependence of time exp[— i(at — mtf>)]. Here, a is, the nor- 
malized in units (M / R 3 ) 1 ^ 2 , oscillation frequency. To lowest 
order for rotating stars, in accordance with the Newtonian 
theory, we expect that the toroidal perturbations of the fluid 
have finite frequencies of order fl (or e). Then from the six 
functions ho, hi, Ho , Hi , H2 , K describ ing the metric pertur- 
bations in the Regge- Wheeler gauge ( jThorne fc Campolat 



taro 1967), only one, ho, is of the same order as the function 
describing the toroidal fluid motions U Q All other pertur- 
bation functions plus the variations of the pressure Sp and 
density Sp are of higher order and thus can be omitted. We 
should point out that this approximation is valid only for the 
study of the r-modes and is not appropriate for the study 
of other fluid or ™-modes. 

In a recent article Lockitch & Friedman (1998) suggest 



that rotation mixes in general all axial and polar perturba- 
tions and introduce the idea of hybrid modes. This mixing 
also eliminates all purely polar modes, but there is still a set 
of purely axial modes. They suggest that this set of purely 
axial modes should not exist for relativistic slowly rotating 
stars at least in the barotropic case. As we show here in the 
specific approximation described above it is possible that 
the r-mode spectrum has a continuous part. 

Using the above assumptions the master equation gov- 
erning quasi- toroidal oscillations is given by (Kojima 199£; 
Kojima 1997|) 



g$ + (zu - fj.) 
where 

and 



«=pr(J-i)(J + 2) 



, 



(9) 



(10) 



(11) 



t The toroidal displacement vector is defined as 

f = (0, u lm sin- 1 e>a , -U tm dg) Y£ 



■ 2 / 

jr w 



167r(p + p)e jr zu = 



(4) 
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q — — i-: {r 4 jzu'^j = 167r(p + p)e x vj 



2m 



(ct — mfi). 



(12) 



(13) 



Equation (^) in the Newtonian limit (j — > 1, g — ► and c<7 — » 
O) reduces to a simple condition for the r-mode frequency 
zu — fi — which is identical with (0) . To this order of f2 the 
eigenfrequency of the r-modes is independent of the radial 
dependence of the eigenfunction <E>(r) (or ho(r)) and in this 
sense the eigenfrequency is infinitely degenerate. 

The master equation (^) is not a regular eigenvalue 
problem since the coefficient (zu — fj.) becomes singular inside 
the star for a certain value of fi. The purpose of this work is 
to study the spectrum of this equation in a rigorous way. 

2.1 Conventions 

The following conventions are used: 

The symbols N, I, C denote the natural numbers (including 
zero), all real numbers greater than zero and the complex 
numbers, respectively. With r we will denote interchange- 
ably some chosen element from 7 or else the identical map- 
ping from 7 to the real numbers. The definition used will be 
clear from the context. For each k € N the symbol C (I, C) 
denotes the linear space of fc-times continuously differen- 
tiable complex-valued functions on 7. 

Throughout the paper Lebesgue integr ation theory in 
the formulation of (Riesz & Sz-Nagy 1955) is used. Com- 
pare with respect to this also Chapter III in (jHirzcbmeTi 
& Scharlau 1971) and Appendix A in (Weidmann 1976). 



Following common usage there is no difference made be- 
tween an almost everywhere (with respect to the chosen 
measure) defined function / and its associated equivalence 
class (consisting of all almost everywhere defined functions 
which differ from / only on a set of measure zero). In 
this sense L 2 C (7, r 4 jj denotes the Hilbert space of complex- 
valued, square integrable functions (with respect to the mea- 
sure r 4 j dr) on the real line. The scalar product < | > on 
Lq (l,r 4 j) is defined by 



<f\9>: 



4 • /■* 

r J f i 



(14) 



for all f,g G L 2 C (/, r 4 j) . Finally, L 2 C (l 2 ) denotes the 
Hilbert space of complex- valued, with respect to the 
Lebesgue measure square integrable functions on the two- 
dimensional interval I 2 . 



3 THE SPECTRUM OF KOJIMA'S MASTER 
EQUATION 

In the following and in the remainder of the paper we con- 
sider the cases / > 1 only. Then, in a first step, we rewrite 
(0) by introducing a new dependent variable ip defined by: 



(15) 



The corresponding representation of $ in terms of tp can 
be performed with the help of special linear independent 
solutions $i,"I > 2 of the differential equation 



77<E> 



('-*)(* + 2) c Ag, 



(16) 



given in the next section. This representation is given by: 



1 

W 



<&2(r) I <&ir' 4 j(pdr' 



$(r) 



where W is the Wronskian of $1 and $2, defined by 



+ $1(0 / $2r' 4 jipdr' 



W := r 4 j ($!$ 2 - $1$ 2 



(17) 



(18) 



Roughly, $1 is an element of Lq (/, r 4 j) for small r and $2 
is an element of L 2 C (/, r 4 j) for large r. Apart from an ir- 
relevant factor such functions are uniquely defined by these 
demands. The use of these functions in the inversion (jl^) of 
( [15] ) is necessary because of the demand (boundary condi- 
tion) that $ is an element of l? c \ I,r A j\. 

Introducing the new variable if into (^) leads to the 
following equation for tp: 



w 



$2(0 



+ (zu(r) - n)tp(r) =0, 



\ $\r tA j<pdr + #i(r) / $ 



/4 • 1 l 

2r j(pdr 



(19) 



for all r £ I. 

In a second step ( Jig| ) is turned into a regular spectral 
problem for the continuous linear operator A on L 2 C (/, r 4 j) 
defined as follows: 

For every / e L 2 C (/, r 4 j) we define Af G L 2 C (/, r 4 j) by 



(Af)(r) 



w 



<J"2(r) / $ir' 4 j/dr' 



+ $i(r) / § 2 r' 4 ifdr' 



+ zo{r)f{r), (20) 



for all r G I. That this indeed defines a continuous linear 
operator on the whole of l? c (l,r 4 f} can be seen as follows. 
First, obviously, since zo is bounded continuous, by 



T^f 



zuf,feL 2 c (l,r 4 j) 



(21) 



there is defined a continuous linear operator on L c (/, r jj 
Second, in the Appendix A it will be shown that from 



(Bf)(r) 



q_ 
W 



$ 2 (r) / $ir' 4 j/dr' 



+ $i(r) / $ 2 r' 4 jfdr' 



(22) 



for all r £ 7 and every / G Lq (7, r 4 j) there is even defined 
a Hilbert- Schmidt operator B on Lq (i , r 4 j) . Hence B is 
not only continuous but in addition compact and Hilbert- 
Schmidt. As a consequence ( poj ) defines a continuous linear 
operator on A on Lq (7,r 4 j) , being equal to the sum of T ro 
and B. 

The determination of the spectrum of Q is now reduced 
to finding the spectrum o~(A) of the continuous linear (non- 
self-adjoint) operator A. Since A is continuous it follows that 
o~(A) is a bounded subset of the complex plane contained in a 
circle around th e origin with the radi us given by the operator 
norm ||A|| of A ( Reed fc Simon 197£ ). The so called essential 
spectrum of the operator is given by the values 
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ZZTq < fJL = 



' + !) 



2m 



(cr - mSl) < n . 



(23) 



(see (Reed fc Simon 1978; Reed fe Simon 1980)) Moreover 



it is known (see e.g. ( Reed fe Simon 1978 ) Corollary 2 of 
Theorem XIII. 14 in Vol. IV) that the essential spectrum t 
is invariant under perturbations by a compact linear oper- 
ator such as B. Hence the essential spectrum a es3 (A) of A, 
which is a part of o(A), is also given by the range of values 
(p3|). The complement a c i isc (A) := a(A) \ a ese (A) (which is 
possibly empty) consists of is olated eigenvalues of finite mul- 
t iplici ty (Reed & Simon 1978). Using an argument of Kojima 



( 1998 ) it follows that these eigenvalues are real (hence cr(A) 
is also real) and are contained in the interval (SI, \\A\\]. 



4 DETERMINATION OF $i AND $2 

Here we distinguish two cases. 

The first uses the following assumptions on p and p : 
Both, p and p are continuous real- valued functions on I sat- 
isfying the Tolman-Oppenheimer-Volkov equations and in 
addition are such that the limits 



lim p(r) and lim p(r) 

r — >0 r— »0 

both exist and 

p(r) = and p(r) = for r > R, 



(24) 



(25) 



are both satisfied, where R denotes the radius of the star. 
[§] Under these conditions the existence of the special solu- 
tions $i,'I > 2 of the differential eq uation ( |l6|) is gi ven by a 
theor e m of Dunkcl (Dunkel 1912) (co mpare also ( Levinson 
19481), ([Bellman 1949Q and (jHille 1969|)) on linear first-order 



systems of differential equations having asymptotic constant 
coefficients at +00. By transformation one gets from this a 
theorem where the singular point is finite. This theorem - 
which for the reader's convenience is given in the Appendix 
A - generalizes well - known results on weakly singular linear 
first-order systems with analytic coefficients. The determi- 
nation of $1 and $2 proceeds now as follows. First from the 
TOV equations it follows that A, j are continuously differ- 
entiable functions on / such that both: 



J — and (26) 

j r 

are continuous as well as Lebesgue integrable near and that 
both 



2 J , 
r — and 



(e A -l) 



(27) 



are continuous as well as Lebesgue integrable near 00. As a 
consequence of the theorem given in the Appendix A follows 
the existence of linearly independent solutions $1, $3 of ( |l6| ) 
satisfyii 



lim r 2 $x( r ) 

r— >0 

lim r l+2 $ 3 (r) 

r—>0 

1+3 



limr i+d $ 3 (r) = -(1 + 2) 



(28) 



and of linearly independent solutions $2, $4 of (|l6|) satisfy- 
ing 

lim r +2 $2(r) = 1 

r — >oo 

lim r i+ V 2 (r) = -(1 + 2) 

r — >oo 

lim r 1_i <I>4(r) = 1 

r — *oo 

lim r 2 - l $' 4 (r) = (29) 

In the next step we conclude that $1, $2 are linear inde- 
pendent. The remaining solutions $3, $4 will be important 
in the proof of the compactness of B (see appendix). 

Up to now everything said in this section is also valid 
for the case 1 = 0. For the following proof of linear indepen- 
dence we have to assume that I > 1. This is assumed in the 
remainder of the paper. 

The proof of the linear independence of $1, $2 proceeds 
indirectly. So we assume on the contrary that there is a non- 
vanishing real a such that $1 = a$2- Using this along with 
(pj[), (p9|), Lebesgue's dominated convergence theorem and 
the monotonous convergence theorem we conclude: 







lim 

n — >oo 



$1 [- (r- 4 ?*;)' + (i - 1)(Z + 2)r 2 je A $j 
-T A j^' x |$, + / r 4 j$fdr 

JR/n 



dr 



+ (l-l)(l + 2) 



R/n 
2 • A=2 , 

r je <P 1 dr 







2 • U!j 

r je <P 1 ar . 



r 4 j<I'i 2 dr + (/-l)(/ + 2) 
'0 ^0 

Where we have made use of the fact that the function 

r 2 je A $ 2 (30) 

is Lebesgue integrable on /. This can be concluded from the 
facts that the function j has a continuous extension onto 
the closed interval [0, 00), that the functions A and j are 
constant for r > R and that 



= 1 
= 1 . 



(31) 



lim r 1 <3?i(r) = 1 

r— >0 



t Defined according to (Reed & Simon 1978) on p. 106 of Vol. IV 



3 These assumptions are satisfied f or instance fo r slowly rotating 
polytropic stellar models (see e.g., (Hartle 1967)). 



lim e A< - r ' 
7 — >o 

lim e A< - r ' 

r — >oo 

All these facts are consequences of the TOV equations and 
the assumptions made on p and p. Hence it follows from 
( pp| ) the contradiction that the function $1 is trivial. As 
a consequence the functions $i,$2 are linear independent. 

The second case considers slowly rotating homogeneous stel- 
lar models. In this case the functions e A and j can b e given 
explicitly in terms of ele mentary functions (see e.g., (CJhan- 
drasekhar & Miller 1974)). It turns out that \,j are contin- 
uous functions and that j is continuously differentiable on 
I \ {R} such that the limits at R of the derivative of j from 
the right and from the left, respectively, differ from each 
other. The functions ( |26| ) are also Lebesgue integrable near 
and ( [27] ) are also Lebesgue integrable near 00. In addition 
( |jl| ) is also valid for this case. As a consequence of the theo- 
rem in the Appendix A follows the existence of continuously 
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different iable functions $1, $2 on I which are two times con- 
tinuously differentiable on I \ \R}, satisfy (|l^) on I \ {R} 
and in addition satisfy (pj|), (]29|). The linear independence 
of these functions can be proved completely analogous to 
the previously considered case. 
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5 DISCUSSION 

In the previous section we showed for a wide class of back- 
ground models for slowly rotating stars that for each ro £ I 
the corresponding /1 = uj(ro) (or rather the corresponding a 
according to (^)) belongs to the spectrum of (^). Further- 
more we achieved a precise definition for the spectrum of ^ 
as the spectrum of the linear operator A. 

This is an interesting new result that needs further 
study. There are still questions to be answered in order to 
understand the effect of the frame dragging induced by gen- 
eral relativity to the spectrum of the rotating stars. For ex- 
ample, one cannot exclu de the possibility that isolated eigen- 
values might also exist (Andersson 1998b). Also, the specific 
form of equation ^ is found under the assumption that the 
r-mode frequencies are of order il. This implies that a few 
terms involving toroidal motions and all the terms related 
to spheroidal motions have been omitted because they were 
either either of higher order or because they contribute via 
higher harmonics (Z± 1, m). This is not the case w ith the hy- 
brid modes found by Lockitch & Friedman (1998) where the 



spheroidal motions have been taken into account. If one in- 
cludes these extra terms the form of the equation will change 
and the effect on the spectrum will be the emergence of the 
hybrid modes, but still the underlying nature of the spec- 
trum can be the same and the techniques applied here can 
be used as well. 

Moreover, in the more general case, there is an imagi- 
nary part for each frequency which corresponds to the damp- 
ing or growth of the fluid motions due to the emission of 
gravitational radiation or due to the viscosity, so it will be 
interesting to examine whether the spectrum will still be 
continuous. 

In other words the discovery of the existence of a con- 
tinuous part of the r-mode spectrum is just the first step 
towards understanding the real nature of the spectrum. Al- 



ready the idea of hybrid modes (Lockitch & Friedmanl998 
Ipser 199^ ) adds new features in the spectrum 



Lindblom 



that have been overlooked in all previous studies. More work 
towards identifying possible isolated eigenvalues of the spec- 
trum should be done and additionally there is need to ex- 
amine if this specific nature of the spectrum is preserved 
when terms of order f2 2 will be included in equation (^|) , the 
coexistence of continuous and discrete part of the spectrum 
is not at all impossible. 

Concluding, we would like to point out that the ex- 
istence of a continuous part of the spectrum might affect 
some of the astrophysical estimations being made for the 
growth time of the r-mode instability and consequ ently all 



the Newtonian estimations being made earlier by ( Anders 
son et al. 1998c 



al. 19981:) 



Kokkotas & Stergioulas 1998 



Andersson et 
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6 H. R. Beyer and K. D. Kokkotas 



APPENDIX A: PROOF OF THE THEOREM 



The variant of the theore m of Dunkel flDun kel 191 4) (com- 
pare also ( [Levinson 1948^ [Bellman 1949| ; pille 196% ) used 
in Section 4 is the following. 

Theorem: Let n £ N; a, to G R with a < to; /1 G JV; 
a M := 1 for = and a M := fi for /i 7^ 0. In addition let 
Ao be a diagonalizable complex n x n matrix and e^, . . . , e' n 
be a basis of C™ consisting of eigenvectors of Ao- Further, 
for each j G {1, • ■ ■ , n} let Xj be the eigenvalue correspond- 
ing to e'j and Pj be the matrix representing the projection 
of C n onto C.e'j with respect to the canonical basis of C n . 
Finally, let Ai be a continuous map from (a, to) into the 
complex n x n matrices M(n x n, C) for which there is a 
number c £ (a, to) such that the restriction of Aijk to [c, to) 
is Lebesgue integrable for each j, k 6 1, n. 



77 (r, r)dr = 



2(2/ + 1) 
(2Z-l)(2Z + 3)' 



k«l 2 



(A7) 



and this expression is integrable over 7 since q is contin- 
uous and has a compact support. Hence the integrability 
of 77 follows from Tonelli's Theorem. Collecting everything, 
we conclude that B defines a Hilbert-Schmidt operator on 
7, 



L c 



,r 4 J) 



Then there is a C 1 map R : (a, to) — > M(n x n, C) with 
Hrrit^oRjk(t) = for each j, fc G l,...,n and such that 
u : (a, to) — * M(n x n, C) defined by 

u(t) := 

' ELi(*o-*r^-(^ + «(*)) --Pj for/i = 
E*=i eajpfAjfto - t)' M ) • (^ + fl(t)) • P, for M / 1 j 

for all t G (a, to) (where 75 is the n x n unit matrix), maps 
into the invertible n x n matrices and satisfies 



(A2) 



for each t G (a, to). 



Now we show that by ( j22| ) for all r G 7 and every / G 
Ijq (7,r 4 j) there is defined a Hilbert-Schmidt operator B 
on Lq (7,r 4 jj. Using the unitary transformation U from 
L 2 c (l,r 4 j) to (7) given by 

Uf:=r 2 y/jf, f£L 2 c (l,r 4 j) (A3) 

it easily seen that this is equivalent to showing that the 
integral operator Int(K') with the kernel function K ' , where 



K'(r,r') 



W 



$2(r)$i(r') for r' < r 



(A4) 



< I > i(r)$2(f') for r' > r. 
for all r G 7 and r' G 7, defines a Hilbert-Schmidt operator 



o n (7). It is well- known (see e.g. (Reed & Simon 1978 



Rccd 



Simon 1980)), that this is equivalent to showing 
uuau i\ is an element of L? c (7 2 ) . Since K' is continuous 
this follows if we can show that \K"^ is integrable over I 2 . 
For this we notice that as a consequence of (pj|), (p9[), there 
are positive real Ci, C2, such that 
(-1 



l*i| 
l*al 



< 
< 



cir 

C2T 



-(1+2) 



(A5) 



From this and since j is (by the TOV equations) bounded, 
the integrability of |7f'| 2 follows as an application of 
Lebesgue's dominated convergence theorem if the integra- 
bility of the following auxiliary function 77 can be shown: 



H(r,r') := \g(r)\ 2 
Now for each r G 7 



r - 2!(r , ) 2 C i+l) for r /< r 

r 2(I+l) (r /j-2I for T , >r 



(A6) 
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